ZEROS OF SYSTEMS OF EXPONENTIAL SUMS AND 
TRIGONOMETRIC POLYNOMIALS 



EVGENIA SOPRUNOVA 

Abstract. Gelfond and Khovanskii found a formula for the sum of the values of 
a Laurent polynomial at the zeros of a system of n Laurent polynomials in (C'*)" 
whose Newton polyhedra have generic mutual positions. An exponential change 
of variables gives a similar formula for exponential sums with rational frequencies. 
We conjecture that this formula holds for exponential sums with real frequencies. 
We give an integral formula which proves the existence-part of the conjectured 
formula not only in the complex situation but also in a very general real setting. 
We also prove the conjectured formula when it gives answer zero, which happens 
in most cases. 



1. Motivation and Summary 

Algebraic geometry is concerned with the study of zero sets of algebraic poly- 
nomials /(x) = X^a"^"^") where a = (ai,...,a„) G Z", x = (xi, . . . , x„) . If 
we allow the exponents a to be real vectors, / becomes a multi-valued function. 
This can be remedied by an exponential change of variables Xj = exp 27rzj , af- 
ter which we obtain a single- valued exponential sum of the form Yla "^o 6xp 2Traz , 
where a = (ai, . . . , a„) G is the vector of frequencies, z = (zi, . . . , Zn) is the 
vector of variables, and az is the standard scalar product. 

Some results from algebraic geometry can be generalized to this wider class of 
functions. For example, Bernstein's theorem states that the number of zeros of a 
generic system of n algebraic equations in (C^)" with a fixed collection of Newton 
polyhedra is equal to n! times the mixed volume of the Newton polyhedra of the 
system. A system of n exponential equations in n variables usually has infinitely 
many isolated zeros, thus one has to study the distribution of these zeros in order 
to obtain finite invariants. O. Gelfond proved that the mean number of complex 
zeros of a system of n exponential sums in n variables with real frequencies whose 
Newton polyhedra have sufficiently general mutual positions (so-called developed 
collection of Newton polyhedra), is equal to n\ times the mixed volume of the 
Newton polyhedra of the system 0. 

O. Gelfond and A. Khovanskii found a formula for the sum of the values of a 
Laurent polynomial at the zeros of a system of n Laurent polynomials in C" with 
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a developed collection of Newton polyhedra [^HH]- This formula splits into two com- 
ponents. One of them is geometrical and reflects the mutual positions of the Newton 
polyhedra of the system, while the other component is expressed explicitly in terms 
of the coefficients of the polynomial and the system. 

An exponential change of variables gives a similar formula for exponential sums 
with rational frequencies. We conjecture that this formula also holds for exponential 
sums with real frequencies. Here is some evidence for this. If the exponential sum that 
we are summing up is identically equal to one, the formula follows from combining 
two results: Gelfond's generalization of Bernstein's theorem [Hj and the new formula 
for the mixed volume J^. In |13j the conjectured formula is proved in dimension 
one. 

The conjectured formula, first of all, implies that the mean value exists. We prove 
the existence not only in the complex situation but also in a very general real setting 
by providing an integral formula for the mean value. 

If the frequencies of the exponential sum that we are summing up are not com- 
mensurate with the frequencies of the system the conjectured formula states that 
the mean value is equal to zero. We show that this is actually true and therefore 
prove the conjectured formula in most cases. For example, if the exponential sum 
that we are summing up is a single exponent explnaz, z G C", a G M", then the 
formula is proved for all values of a except for a countable set in M" . 

Our arguments represent a combination of real analytic geometry and ergodic 
theory. They are based on two theorems from completely different parts of mathe- 
matics: the cell-decomposition theorem for subanalytic sets (Appendix B) and Weyl's 
equidistribution law for multidimensional trajectories in the real torus (Appendix 
A). Weyl's equidistribution law for one-dimensional trajectories is a classical the- 
orem which was published in 1916 (see ^01 )• The proof of the corresponding law 
for multidimensional trajectories in jT2] is a direct generalization of Weyl's original 
argument. 

Remarks. Similar ideas first appeared in Weyl's papers ^1] and JS| where he 
solves the mean motion problem. S. Gusein-Zade and A. Esterov were dealing with 
related questions and were using close techniques in ^ ll()| |3] to prove the existence 
of the mean Euler characteristic and mean Betti numbers of level sets and sets of 
smaller values of a quasiperiodic function. Their motivation is completely different 
from ours: it comes from the analysis of some models of chaotic behavior appearing 
in quasicrystal structures. See also a recent paper jlj of V. I. Arnold for a close 
discussion. 

Acknowledgements. This paper is a part of the author's Ph.D. thesis defended 
in spring 2002. I would like to thank my thesis advisor Askold Khovanskii for stating 
the problem and for his constant attention to this work. 
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2. The Gelfond-Khovanskii formula and its conjectured 
generalization to the case of exponential sums with real 

frequencies 

Let Ai, . . . , A„ be convex polyhedra in M" and A = Ai + • • ■ + A„ be their 
Minkowski sum. 

Definition 2.1. A collection of faces {Tj : Tj C Aj,j = l...n} is coordinated if 
there exists a nonzero linear function on R" whose maximum on the polyhedron A-,- 
is attained exactly on the face Fj C Aj for all j = 1 . . . n . A collection of polyhedra 
Ai, . . . , An is developed if none of the polyhedra is a vertex, and in any coordinated 
collection of faces there is at least one vertex. 

Fix a system 

Pi(z) = ••• =P„(z) =0, zGC" 
of n Laurent equations with a developed collection of Newton polyhedra Ai, . . . , A„ . 
The Gelfond-Khovanskii formula states that the sum of the values of a Laurent 
polynomial Q over the zeros of the system in (C^)" is equal to 

a 

where the summation is performed over the vertices a of A = Ai + • • • + A„ , ka 
is the combinatorial coefficient that corresponds to the vertex a (the combinatorial 
coefficients are integers that reflect the mutual position of the Newton polyhedra 
Ai, . . . , An of the system, see (SHE] for the definition), and Ca is an explicit Laurent 
polynomial in the coefficients of P = Pi • • • P„ and Q . This result was announced 
in [S] , and a proof was given in [HI . 

If we allow the exponents of the Laurent polynomials to be real vectors, we ob- 
tain multi- valued functions. After an exponential change of variables they become 
single- valued exponential sums with real frequencies. Now we make some prepara- 
tions before we formulate the conjectured generalization of the Gelfond-Khovanskii 
formula to the case of exponential sums. 

Let A be a finite set in M" . An exponential sum with the spectrum A is a function 
P : ^ C of the form 

F{z) = Cq, exp 2-Kaz, 

aeA 

where the summation is performed over the frequencies a G A, z = {zi, . . . , Zn) G 
C", Cq are nonzero complex numbers, and az is the standard scalar product. The 
Newton polyhedron of an exponential sum is the convex hull A(A) of its spectrum A. 
Fix a system of exponential sums 

(2.1) Pi(z) = ... = p„(z) =0, zGC" 

with a developed collection of Newton polyhedra Ai,...,A„ in M". Due to the 
assumption that the collection of the polyhedra is developed, there exists P > 
such that all the zeros of the system lie in a strip Sr x Im C"' , where Sr C Re C" is 
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a ball of radius R centered at the origin. This implies that all the zeros are isolated 
(see Theorem 15. ip . 

Let G be another exponential sum with real frequencies. The sum of the values at 
the zeros of a system in the exponential case is replaced with the result of averaging 
G over the zeros of the system along the imaginary subspace. Let C Ln C" be a 
bounded set with nonzero volume. For A € M, define Sq{X) to be the sum of the 
values of G at the zeros of (|2.1|) (counting multiplicities) that belong to the strip 
M" X AO C C" . 

Definition 2.2. The mean value of G over the zeros of the system (|2.1j) is the 
limit of /^^(A)/ Vol(Ar2) , as A approaches infinity. 

The Minkowski sum A = Ai + • • • + A„ is the Newton polyhedron of the product 
F = Fi - ■ ■ Fn- The exponent exp 2Traz that corresponds to a vertex a in A appears 
in F with a nonzero coefficient da ■ Let F = F/ {da exp 2'Kaz) . The constant term 
of the exponential sum F is equal to one. We define the exponential series for 1/F 
by the formula 

l/F = l + {l-F) + {l-Ff + ■■■ . 

Since each exponent appears with a nonzero coefficient in a finite number of terms, 
the coefficients of this series are well-defined. Let Ca be the constant term in the 
formal product of this series and {\/da )exp(-27raz)Gdet(|f). 

Conjecture 2.3. The mean value of an exponential sum G over the zeros of 
the system (|2.1|) is equal to 

(_27r)" 

where the summation is performed over the vertices a of the Minkowski sum A = 
Ai + • • • + A„ , Ca are the constant terms of the series as defined above, and ka are 
the combinatorial coefficients. 

As we mentioned in the summary, this formula is proved in the following three 
cases: 

1. The frequencies a are rational (the formula is obtained by an exponential 
change of variables from the algebraic case). 

2. G = 1 (the formula follows from combining two results: Gelfond's gener- 
alization of Bernstein's theorem jSj and the new formula for mixed volume 

El)- 

3. n = 1 (see [13]). 

In the exponential case even the existence of the mean value is not obvious. 
We give an integral formula for the mean value which proves the existence not 
only in the complex situation (Theorem 15. 2|) but also in a very general real setting 
(Theorem 14. 3j) . 

Assume that the frequencies of the exponential sum that we are summing up are 
not commensurate with the frequencies of the system. Then there is no constant term 
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in the series defined above, that is, all Cq, are equal to zero, and the conjectured 
formula states that the mean value is equal to zero. Using the integral representation, 
we show in Corollary that this is actually true. This proves the conjectured 
formula in most cases. For example, if the exponential sum that we are summing up 
is a single exponent exp I'Kaz, then we proved the formula for all values of a except 
for a countable set in R" . 

3. Averaging over the isolated intersections of a subanalytic set 

WITH A dense orbit ON 

Here we present a construction and formulate a theorem on which all our results 
are based. Let y be a subanalytic subset of the real torus and O a dense 
multidimensional orbit on . We prove an integral formula for the mean value of a 
bounded subanalytic function on V over the isolated intersections of the set V and 
the orbit O . This formula, in particular, implies that the mean value always exists. 

Let X = (xi, . . . , Xn) G and (/J = (991, . . . , (/Jtv) G = M^/Z^ . A linear map 
$ : ^ defines an action of on by 

X : 93 1-^ (93 + <I>(x)) modZ^. 

Let O be an orbit of this action. We require that the orbits are dense in the torus, 
which means that there are no integral vectors orthogonal to the n-plane $(R"') C 
(Appendix A). We will also assume that <I> is injective, that is, the orbits are 
n-dimensional. 

Let y be a subanalytic subset of and T a bounded subanalytic function on 
V. Set V := <^-^{V n O) C M", and define a function T(x) = f{<^{x)) on V. Let 
O be a bounded subset of with nonzero volume. For A > 0, define ^^(A) to be 
the sum of the values of T at the isolated points of V that belong to . 

Definition 3.1. The mean value Mq of T over the isolated intersections of V and 
O (in the topology of O) is the limit: 

A^oo Vol(Ar2) 

Let p : ^ L be the orthogonal projection to the linear subspace L C 
orthogonal to $(M"). Define MAr_„(y) to be the smooth (A^ — n) -dimensional part 
of V , transversal to the orbit. Here is the precise three step definition (see Figure 1 
below where we assume that the direction of the orbit is vertical): 

(1) Let Vo be the set of all points ip G V such that the intersection of V with 
the n-plane through ip parallel to the orbit is locally just the point 97 itself. 

(2) Let Vi be the set of all points (p G Vq such that Vq is a C"*^ -manifold of 
dimension N — n in a small neighborhood of if. 

(3) Let M]y-n(y) be the set of all points ip £ Vi such that the projection p is 
regular at each of these points. 
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(1) (2) (3) 

Figure 1. Definition of MAr_„(l/) 



Let = ^(ej) , where ej is the j-th vector in the standard basis for M". The 
quotient of the standard volume form on R-^ by the image under $ of the stan- 
dard volume form on R" defines a volume form uj = dip{. . . , , . . . ,^") and an 
orientation on L. Then M^-niV) is a -manifold, and the projection p defines 
an (iV — n)-form p*uj on Mjv-n(^) • 

Definition 3.2. The O -transversal T -weighted volume of V is the integral of the 
form Tp*iJ over Mjv_„(y). 

Definition 3.3. We say that the fractal dimension of a set X is less than or equal 
to n if there exist constants c and eo > such that for any < e < eg , the set 
X can be covered by no more than ce~" balls of radius e. We say that the fractal 
dimension is less than n if it is less than or equal to n — 6 for some positive 6 . 

Our main result is the following theorem: 

Theorem 3.4. Let O he an n -dimensional dense orbit in the torus = M^/Z^, 
n < N , V a subanalytic subset of , and T a bounded subanalytic function on V . 
Let 0, be a bounded subset of M" with nonzero volume whose boundary has fractal 
dimension less than n . ^ ^ 

Then the mean value Mq of T over the isolated intersections of V and O always 
exists and equals the O -transversal T -weighted volume of V . 

We prove this theorem in Section 6. 

4. Zeros of systems of trigonometric polynomials 

We apply Theorem 13.41 to prove that the mean value of a quasiperiodic trigono- 
metric polynomial T over the isolated points of a quasiperiodic semitrigonometric 
set V (a set described by quasiperiodic trigonometric polynomials, see the precise 
definition below) always exists. We obtain an integral representation for the mean 
value from which we conclude that in the case when the frequencies of T are not 
commensurate with the frequencies of the quasiperiodic trigonometric polynomials 
that describe V , the mean value is equal to zero. 
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4.1. Definitions. Let x € M". A quasiperiodic trigonometric polynomial is a 
trigonometric polynomial of the form 

p 

T{x) = Cfc cos InakX + dk sin InokX 

k=l 

with real coefficients Ck,dk € M and real frequencies G M". Here a^x is the 
standard scalar product. 

A subset of M" is called quasiperiodic semitrigonometric if it is described by 
quasiperiodic trigonometric equations and inequalities. More precisely, we say that 
V C M" is quasiperiodic semitrigonometric if it can be represented in the form: 

1^ = u n 

i=ij=i 

where each Vij is either {x G |Tjj(x) = 0} or {x G M" \Tij{x) > 0}, for some 
quasiperiodic trigonometric polynomials Tjj . 

Let T be a quasiperiodic trigonometric polynomial, V a quasiperiodic semitrigono- 
metric set in described by Tij , and Q a bounded subset of M" with nonzero 
volume. For A > 0, define ^^(A) to be the sum of the values of T at the isolated 
points of V that belong to AO. The mean value of T over the isolated points 
of V is defined by 

Mn = hm 

A^oo Voi(Ar2) 

4.2. Averaging a quasiperiodic trigonometric polynomial over the isolated 
points of a quasiperiodic semitrigonometric set. Consider the set of all fre- 
quencies of T and Tij . This set generates a subgroup A of (M" , +) , which is a direct 
sum of infinite cyclic subgroups: 

A = (^i) e • • • e (Aat). 

That is, there exist Ai, . . . , Aj^ G with no integral relations, such that each of 
the frequencies is an integral combination of ^i, . . . , ^4^^ • 

Define a linear map <I> : M" ^ by <I>(x) = {Aix, . . . ,Anx) . Then R" acts on 
= rN/Z^ by 

X : (fh^ + <I>(x)) modZ^. 

Let O be the orbit of this action through the origin. Since each of the frequencies 
is an integral combination of ^i, . . . , A^v , the functions T and Tij are the restrictions 
to O of some trigonometric polynomials T and Tij defined on . For example, 
cos(27rQ!x) is the restriction of cos 27r ^ mj^jj , where a = "^ruiAi for rrii^'L. 

To apply Theorem 18.41 here we need the orbit to be dense (Lemma 14.11) . the 
mapping $ to be injective (Lemma 14. 2|) . and the dimension n of the orbit to be 
less than (if $ is injective then n < N ; N = n corresponds to the periodic case 
which we treat in R.emark 14. 4|) . 
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Lemma 4.1. There are no nonzero integral vectors orthogonal to the plane of the 
orbit O . Therefore, O is dense in the torus . 

Proof. Let A be the N xn matrix whose rows are the generators Ai, . . . , of A. 
We denote the columns of this matrix by A^, . . . , A"" . These columns generate the 
n-plane <I>(R") of the orbit O . If there exists a nonzero integral vector (ki, . . . , k^) 
orthogonal to O then it is orthogonal to the vectors A^, . . . , A^ . This implies that 
kiAi + - ■ ■ + kNAN = 0, i.e. there exists a nontrivial linear combination with integral 
coefficients on the generators Ai, . . . , Aj\[ of A, which gives a contradiction. By 
Appendix A the orbit O is dense in the torus . □ 

Lemma 4.2. If the dimension of the orbit O is less than n, then the set V described 
by the Tij has no isolated points, and the mean value Mq is equal to zero. 

Proof. If the dimension of the orbit is less than n, then A^,...,A^ are linearly 
dependent, and the solution space of the system A^xi + • • ■ + A^Xn = is nontrivial. 
If Tij(x) = 0, then Tij{x + y) = for each y in the solution space. This means that 

V has no isolated points, and the mean value is equal to zero. □ 

If the orbit is ndimensional and n < N , then M" is mapped bijectively to a 
dense n -dimensional orbit on the real torus and trigonometric polynomials Tij 
regarded as functions T,j on define a set V . The intersection of the set V and 
the orbit O is the semitrigonometric set V . We are adding up the values of T over 
the isolated (in the topology of the orbit) points of intersection of the set V and 
the orbit. Here is the main result of this chapter which now follows directly from 
Theorem 13.41 

Theorem 4.3. Let V be a quasiperiodic semitrigonometric subset of M", T a 
quasiperiodic trigonometric polynomial. Let V , T , O be the corresponding set, func- 
tion on , and orbit, as constructed above. Let ^ be a bounded subset of M" with 
nonzero volume, whose boundary has fractal dimension less than n . 

Then the mean value always exists and is equal to the O -transversal 

T -weighted volume of V , assuming that the dimension of the orbit is n, and n < N . 
If the dimension of the orbit is less than n , the mean value is equal to zero. 

Remark 4.4. If the orbit has the same dimension as the torus, i.e., N = n, the set 

V C and the function T have n linearly independent periods A^, . . . , A^ , which 
define a torus W/{ZA^ • • • ZA^) . The mean value Mq in this case is equal to 
the sum of the values of T at the isolated points of V in this torus, divided by the 
volume of the torus. 

Remark 4.5. This theorem holds and our proof works if y is a quasiperiodic 
subanalytic subset of M", that is, V is the intersection of a subanalytic subset V of 
with some dense n -dimensional orbit, and T is the restriction to this orbit of 
some bounded subanalytic function on V . 
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Corollary 4.6. // all the multiples ka, /c S Z \ {0} of each frequency a of T do 
not belong to the subgroup A of (IR"',+) generated by the frequencies of Tij , then 
the mean value Mq is equal to zero. 

Proof. It is enough to check this statement for T with a single frequency a, i.e., 
for T[x) = ccos27rQx + dsinlTiax . The set of all frequencies of the trigonometric 
polynomials Tij describing the set V generates a group A! in which is the 

direct sum of infinite cyclic subgroups: 

M = {Ai) e • • • © (^7V-i). 

Since no multiple of the frequency a of T belongs to A! , the group A obtained 
from A' by throwing in a, is the direct sum: 

A = {Ai) • • • © {An-i) © (a). 

Therefore, raising the situation to the torus M^/Z^ , we obtain a function 
T{(Pn) = ccosifN + ds'mipN that depends only on the variable (pN , and a set 
V , which is described by the functions Tij that do not depend on ipN ■ 

By Theorem l4.31 the mean value is equal to the O -transversal T-weighted volume 
of V. Since the Tij do not depend on ipj^ , the manifold M^_n{V) (which is the 
smooth N — n dimensional, transversal to the orbit O, part of V) is the direct 
product V X , where V is a manifold in M^-VZ^~i , and = R/Z. The form 
uj = p*dipi . . . dipN{- . . , . . . , A"-) can be written as uji A dip^ where uJi is a form 
on V' . Therefore, the mean value is equal to 

/ _ Too = ( / {ccos2Trip]\f + dsm2'KipN)dip]\[ \ LJi = 0. 

JMN-niV) JV \JS^ J 

□ 

5. Zeros of systems of exponential sums 

In this section we apply our main result (Theorem 13. 4() to the complex case, that 
is, to computing the mean value of an exponential sum over the zeros of a system 
of n exponential sums in C". We show that this mean value always exists. From 
the integral representation for the mean value that we obtain, we deduce a proof 
of Conjecture 12.31 in the case when the frequencies of the exponential sum are not 
commensurate with the frequencies of the system. 

5.1. Definitions. Fix a system 

(5.1) Fi(z) = --- = F,(z) = 0, zGC" 

of exponential sums with a developed collection of Newton polyhedra. Let T be 
a family of systems of n exponential sums with the same collection of spectra as 
(|5.H) . and whose coefficients have the same absolute values as the corresponding 
coefficients in (|5.1j) . The following theorem is a particular case of a result proved by 
Gelfond 0E|. 
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Theorem 5.1. There exists R > such that all the zeros of the systems from the 
family T belong to the strip Sr x Im , where Sr C Re is a hall of radius R 
centered at the origin. This implies that all the zeros of the systems from T are 
isolated. 

Let G be an exponential sum with real frequencies, and 0, a bounded subset of 
ImC" with nonzero volume. For A > 0, let Sn{X) be the sum of the values of G at 
the zeros of (|5.1() (counting multiplicities) that belong to the strip x XQ c C". 
Let Sfi{X)-^ be the sum of the values of G at the zeros (not counting multiplicities) 
of the system (|5.1() of multiplicity at least k that belong to the strip M" x AO C C" . 
By Theorem 15.11 Sq{X) and S'q(A)-'^ are well-defined. 

Define 

A^oo Vol(Ail) A^oo Vol(Ar2) 

Theorem 15.21 below provides an integral formula for computing the mean values 
Mq , similar to the formula in the real case (Theorem 14. In particular, The- 

orem 1^21 implies that the mean value Mq always exists. Another consequence of this 
integral formula is that in the case when the frequencies of G are not commensurate 
with the frequencies of the exponential polynomials from the system, the mean value 
Mq is equal to zero fCorollarv I5.,'-{|) . 

5.2. Averaging an exponential sum over the zeros of a system of exponen- 
tial sums. As in the real case, we deduce an integral formula for the mean value 
from Theorem K-i.41 The exponential sums Fj , G are quasiperiodic along the imag- 
inary subspace of C" . Since the system (|5.1j) has a developed collection of Newton 
polyhedra, the real parts of the zeros of this system belong to some ball Sr C Re C" . 
Therefore, we can restrict our attention to this ball, and think of the Fj and G as 
of functions quasiperiodic along the imaginary subspace and periodic along the real 
subspace. This allows us to consider C" as an orbit on some real torus and apply 
Theorem 13.41 to this situation. 

As before, the set of all frequencies a of the exponential sums Fj , G generates a 
subgroup A in (M",-!-), which is a direct sum of infinite cyclic subgroups: 

A = (ai) © • • • e (oat). 

In other words, there exist ai,. . . , a^v £ 1^" with no integral relations, such that 
each of the frequencies a G ^ is an integral combination of ai, . . . ,a]\f ■ 
Define a linear map 

$i:ImC"^]R^ by $i(y) = (aiy, . . . , ajvy), y SlmC". 

According to Theorem 15.11 the real parts of the zeros of the system (|5.1|) belong 
to the open ball Sr C ReC" centered at the origin. Let $o : ReC" ^ R" be a 
composition of a shift and a rescaling that maps Sr inside the unit cube. We now 
think of Fj and G as of functions periodic along the real subspace. 
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Then $ = $o x ^i maps C" = Re C" x Im C" to R" x . The map ^> defines 
an action of C" on T"+^ by 

(99 + $(x,y))modZ"+^, G T"+^, xeReC", y G Im C". 

Let O be the orbit of this action through the origin. Since each of the frequencies a 
is an integral combination of ai, . . . , a„ , the sums Fj , G are the restrictions to O 
of some functions Fj , G defined on T^^^ . We note that these functions may only 
be non-analytic at the points (p where at least one of the coordinates , . . . , is 
equal to zero. 

The set of all zeros of the system ()5.1() of multiplicity at least k is an analytic set 
defined by exponential equations Fj{z) = and a few more exponential equations 
whose frequencies belong to the subgroup A of (see, for example, j^). 

Let V^^ be the set of aU zeros of the system of multiplicity greater than or 
equal to k . Consider the system of exponential sums that defines V-^ . Regarded 
as functions on Tr"'^'^ , these sums define a subset V-^ of the torus T""''^ . We now 
state the main result of this section. 

Theorem 5.2. Consider a system of exponential sums with real frequencies 
(5.2) Fi(z) = ■•• = F„(z) =0, zgC" 

with a developed collection of Newton polyhedra. Let G he an exponential sum with 
real frequencies, and Q a hounded subset of M" with nonzero volume whose boundary 
has fractal dimension less than n. Let V-^, G , O be the corresponding subsets of 
T""*"^, function on ']r"+^^ and orbit, as constructed above. 

Then if n < N the mean value is equal to the G -weighted O -transversal 

volume of V-^ . The mean value Mq is equal to the sum + • • • + for 

some k. 

Proof. We first show that the dimension of the orbit O is 2n. It is enough to 
show that the dimension of the plane $i(Im C") is n. This plane is generated by 
the columns of the matrix of the linear map . The rows of this matrix are the 
generators ai, . . . , oat of the group A. Thus the dimension of the plane $i(Im C") 
is equal to the dimension of the linear space generated by the frequencies ai , . . . , oat . 
Since the collection of the Newton polyhedra of the system (|5.2|) is developed, its 
frequencies generate . We conclude that the dimension of the orbit is n . 

Similarly to the real case, it is easy to see that there are no integral vectors 
orthogonal to the orbit O . Therefore, the orbit O is dense in the torus Y^^^ (see 
Lemma 14. 

To derive Theorem 15. 21 from Theorem 13 . 41 we show that V-'^ is an analytic subset 
of T"^^ and G is analytic on V-'^ . Each of V-^ is defined by a system of equa- 
tions analytic everywhere on T'^+^ but possibly at the points where at least one 
of the coordinates ipi, . . . ,ipn is equal to zero. Notice that V-'' is contained in the 
"cylinder" ^q{Sr) x T^. Indeed, the intersections of V with the orbit through the 
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origin belong to this cylinder since the zeros of the initial system lie in the strip 
Sr X Im C" . Restricting the functions Fj to a shifted orbit we obtain functions 
that differ from Fj by a shift along the imaginary subspace. By Theorem 15.11 the 
zeros of this new system belong to the strip Sn x Im C", and we conclude that 
\/>fc (2 ^i(^Sr) X T^. Therefore, V-^ is locally defined by a system of analytic 
equations, and G is analytic on V-^. 

To show that Mq can be computed by the formula Mq = + • • • + 

for some k, we need to prove that the sets V-'' are empty, starting from some k. 
Assume that there exists a sequence = + ii/k, € Sr, such that belongs 
to V-^. Raising this sequence to Sr x T^, we obtain a sequence z^ in Sr x 
such that ifc belongs to V-''. Since Sr x is compact, there exists a limit point 
zq & Sr X . This point zq belongs to each V-'^. Although this point does not 
necessarily belong to the orbit Sr x <I>i(Im C^/Z^) , it belongs to some shifted orbit. 
The restrictions of the functions Fj to this shifted orbit define an exponential system 
whose zeros are isolated and belong to Sr x Im C" fTheorem l5.1j) . Therefore, zq is 
a zero of the shifted system of infinite multiplicity, which is impossible since all the 
zeros of this system are isolated. We proved that all V-^ are empty starting from 
some k. 

Now the theorem follows from Theorem 13.41 □ 

The proof of the following corollary is a repetition of the proof of Corollary 14.61 

Corollary 5.3. If all the multiples ka, k & 'L \ {0} of each frequency a of F do 
not belong to the subgroup A of (M",+) generated by the frequencies of the Fj , then 
the mean value Mq, is equal to zero. 

6. Proof of Theorem 13.41 

Proof Recall that O = {a + ^>(M'')) mod , where a G T^, and $ : ^ 
is a linear injective map. We send the torus to the cube / = [0, 1)^ C by 
cutting along (pj = 1. Choose a new basis {B\ . . . , S^"", A^,...,A'^} in 
whose last n coordinate vectors are the images of the standard basis in R" under 
$, A^ = ^(cj) . Note that the vectors A^ , . . . , A" generate the plane of the orbit 
$(R"') . Let the remaining vectors B^, . . . , B^~'"' form a basis of the (A^ — n)-plane L 
orthogonal to the orbit, so that the determinant of the change of coordinates matrix 
from the standard basis in R-^ to this new basis is equal to one. 

Let X = {xi, . . . , Xn) be the coordinates along the orbit, and y = (yi, . . . , yN-n) 
the coordinates in the orthogonal plane L. We choose a cube Q = Qy x whose 
edges go along the new coordinate axes and I <Z Q . We now regard ^ as a subset 
of / C Q and T{ip) as a function T(y, x) on V <Z F 

The set y is a subset of Q , subanalytic in RP^ , and the graph of the function 
T is a subset of Q x R subanalytic in RP^"*"^ . Therefore, by the cell decomposition 
theorem (Theorem 18. 3j) . there exists a decomposition of R^ with the fixed order of 
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coordinates (y, x) into finitely many cells such that ^ is a union of cells and T is 
on each of these cells. 

Since the cell decomposition respects the standard projection, each cell is either a 
part of the closure of another cell, or is disjoint from it. Only (ii, . . . , iN-n, 0, . . . , 0)- 
cells that do not belong to the closure of a cell of higher dimension may have isolated 
intersections with the orbit. We will show later that cells of dimension less than 

— n do not contribute to the mean value. Therefore, we can restrict our attention 
to (1, . . . , 1, 0, . . . , 0)-cells with N — n ones and n zeros that do not belong to a 
closure of a cell of higher dimension. 

Let C be one of such cells. Then 



where H = {Hi, . . . , Hn) for some functions Hi, . . . , Hn , and D is a (1, . . . , 1)- 
cell in Qy . 

Now we define a counting function g on the torus . Through each point a C 
we draw an n-plane parallel to the orbit. If tp belongs to (a + <I>(i3r))/Z^ , where 
Br C is the ball of radius r centered at the origin, we define g{'p) = T[a)/vr, 
where Vr is the volume of the ball B^ C M" . (Here r is chosen so small that 
the mapping x i— > ^>(x)/Z^ restricted to Br is one-to-one.) Otherwise we define 



Figure El illustrates the definition of g in the case when the cell C is one-dimensional. 
Through each a € C draw a vertical segment of length 2r centered at a. This 
segment is the ball a + Br from the general construction. Define g to be identically 
equal to T{a)/2r on that segment. 



C = {{y,x) \ yGD, x = H{y)}, 




T{y,H{y))/vr, iiyeD, xeBr 



0, otherwise. 



X 




Figure 2. Definition of the counting function g 
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Since the function T is continuous and bounded on C and the functions defining 
C are continuously differentiable, g is Riemann integrable on the torus T-^. We 
apply Weyl's equidistribution law for multidimensional trajectories (Theorem [7^: 

(6.1) lim \ / gmx))dx = [ g{^)dip. 

A^oo Vol(Ai2) Jxn JjN 

Let us show first that the left hand side part of this identity is exactly the contri- 
bution of the cell C to the mean value Mq . a + Bj. , sl ball centered at the point 
a € M" , lies entirely in then by the definition of g 

g{<!>{x))dx = T{a). 

a+Br 

It remains to take care of the points a such that a + Br does not lie entirely 
in Ari. The dimension of the boundary of is less than or equal to n — 6 for 
some 6 > 0. Hence the number of balls of radius r that are required to cover the 
r-neighborhood of the boundary of AO (i.e. the set of all points whose distance to 
the boundary is less than r) is bounded by cX^~^ for some constant c. The number 
of isolated points of V in each of such balls is bounded by the number of cells in 
our cell decomposition. Therefore, the result of adding up T over the isolated points 
of V , whose distance to the boundary of AO is less than r, is bounded by cqA""^ 
for some constant cq. When we divide the result by the volume of AO, which is 
proportional to A" , the quotient approaches as A approaches infinity. Therefore, 
we can disregard the points that are close to the boundary of AO while computing 
the mean value Mq . We have proved that 



A^oo Vol(AO) Jxn 



is the contribution of the cell C to the mean value Mq . 

We now work on the rewriting of the right hand side part of (|6.H1 . 



g{(f)d(p = / g{y,x)dydx= / / g{y,x)dx\dy 
T^v Jq jQy \Jq^ / 

fiy,Hiy))dy= [ fiip)p*d^i...,A\...,A^)= [ f{^)p*uj. 

ID JC Jc 

Next we note that cells of dimension less than — n do not contribute to the 
mean value. Indeed, for each of such cells we construct a counting function g and 
repeat the above argument. It follows that the contribution of each of such cells is 

g{ip)dip, 

JN 

which is equal to zero since g is zero on a set of full measure in . 

The union of all (ii, . . . ,iAr_„,0, . . . ,0)-cells that do not belong to the closure 
of a cell of higher dimension is the {N — n) -dimensional manifold M^_„(y), up 
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to possibly a few cells of smaller dimension. Those cells of smaller dimension do 
not contribute to the mean value. We conclude that the mean value is equal to the 
T-weighted O -transversal volume of V . □ 

7. Appendix A: Weyl's equidistribution law for dense 

MULTIDIMENSIONAL ORBITS ON THE REAL TORUS. 

Classical Weyl's equidistribution law JHl states that the time-average of a Rie- 
mann integrable function along a one-dimensional dense orbit on a real torus coin- 
cides with the space average. Here we state a version of Weyl's law for multidimen- 
sional orbits. 

Let X G R" , G = M^/Z^ . A hnear map $ : M" ^ defines an action of 
M" on T'^ by ip i—f {ip + <I>(x)) modZ^ . Let O be an orbit of this action. 

Theorem 7.1. Let f be a Riemann integrable function on the torus , O an 
orbit on through an arbitrary point ipo € T-^ such that there are no nonzero 
integral vectors orthogonal to O , and a subset of with nonzero volume. 
Then 

.r i/\n\ / f{^o + ^{x))dx^ j f{<f)d(f, as A ^ oo. 
Vol(AO) Jxn JrN 

The proof of this theorem in {12) is a direct generalization of Weyl's original 
argument for one-dimensional orbits. 

Remark 7.2. This theorem implies that an orbit is dense if and only if there are 
no nonzero integral vectors orthogonal to the plane of the orbit. 

8. Appendix B: Cell-decomposition theorem 

The cell-decomposition theorem for subanalytic sets states that a subanalytic set 
can be partitioned into finitely many cells, which are subanalytic subsets of especially 
simple form. We will follow the formulation of the cell-decomposition theorem as 
presented in [^J. 

Definition 8.1. Let {ii, . . . ,in) be a sequence of zeros and ones of length n. An 
{ii, . . . , in)-cell is a subset of M" obtained by induction on n as follows: 

(i) a (O)-cell is a point in M, a (l)-cell is an interval (a, 6) € M; 

(ii) suppose (ii, i„) -cells are already defined; then an (zi, . . . , i„, 0)-cell is 
the graph r(/) of a function / on an (zi, . . . , i„)-cell provided that 
r(/) C M"^"*^ is subanalytic in MP"+^ ; further, an (ii, . . . , i„, l)-cell is a set 

:= {(.x,r) eX xR: f{x) <r< g{x)}, 

where X is cin ('^i, . . . ,i7T,)-CGll, functions fiQ circ functions on JC wliosG 
graphs T{f),T{g) C M"+^ are subanalytic in , and / < g on X. The 

constant functions f{x) = — oo, g{x) = +oo are also allowed. 
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Definition 8.2. A decomposition of M" is a special kind of partition of M" into 
finitely many cells. The definition is by induction on n: 

(i) a decomposition of M is a collection 

{(-oo, ai), (ai, 02), . . . , (afc, +00), {ai}, . . . , {0^}} 

where ai < ■ ■ ■ < Ok are points in R. 

(ii) a decomposition of M""*"^ is a finite partition into cells such that the set of 
projections of these cells to the first n coordinates forms a decomposition of 
W. 

Theorem 8.3. (The cell decomposition theorem) Let X be a subset of M" 
subanalytic in MP"- . Let f be a function on X whose graph in M"+-^ is a subanalytic 
subset of MP"^^ . Then there exists a (finite) decomposition V of M" such that X 
is a union of cells in V and the restriction flB-B^M. to each cell B ^ V is . 
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